SHOCK  SENSITIVITY,  A.  ?*0P£3TY  OF 

many  &ma$ 


ATTENTION 


NOLTP.  70-73 


SHOCK  SENSITIVITY,  A  FROPERTY  CP  MANY  ASPECTS 


Donna  Price 


ABSTRACTS  (U)  Shock  sensitivity  of  an  explosive  includes  its 
threshold  for  propagation  of  steady  state  detonation  at  the  critical 
diameter  (d  )  as  well  as  its  numerous  thresholds  for  initiation  of 
detonation  under  different  transient  conditions.  Data  for  TNT 
charges  of  differing  degrees  of  homogeneity  are  used  to  show  a 
continuous  variation  of  dQ  with  the  Initiating  pressure  P^  measured 
with  the  NOL  large  scale  gap  test.  An  example  of  a  critical  curve 
is  constructed  in  the  pressure-time  plane;  it  runs  from  the  threshold 
conditions  at  dQ  through  those  at  P^.  It  Is  suggested  that  all 
other  initiating  pressures  measured  for  the  same  charge  would  also 
fall  on  this  curve  which  illustrates  the  relationship  between  the 
different  threshold  values. 
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SHOCK  SENSITIVITY,  A  PROPERTY  OF  MANY  ASPECTS 

The  "shock  sensitivity"  of  an  explosive  cannot  be  described  by 
a  single  valued  parameter.  Shock  sensitivity  encompasses  all  of  the 
complex  reactions  of  the  explosive  to  many  different  shock  conditions. 

•  Such  reactions  are  manifested  by  the  phenomenon  of  failure  diameter 

for  propagation  of  detonation  as  well  as  by  the  critical  initiation 
,  pressures  for  shock  to  detonation  transition  observed  in  the  numerous 

gap,  booster,  projectile,  and  wedge  experiments.  A  complete  descrip¬ 
tion  of  an  explosive's  sensitivity  to  shock  should  encompass  there¬ 
fore,  not  only  the  thresholds  for  initiation  under  transient  conditions 
but  also  the  threshold  for  propagation  under  steady  state  conditions. 
Moreover,  it  should  account  for  the  large  effect  of  the  physical 
state  of  the  explosive  on  both  thresholds. 

The  relationships  between  these  many  aspects  of  shock  sensitivity 
are  complex  and  cannot  be  simply  stated.  This  situation  has  led  to 
a  general  confusion  about  the  meaning  of  different  measurements.  It 
is  the  purpose  of  this  papar  to  clarify  the  situation  by  describing 
the  relationships  between  the  conditions  at  the  critical  diameter 
and  those  at  the  critical  initiating  pressures  (gap  test)  and  then 
by  using  measured  values  for  one  explosive  in  a  number  of  different 
physical  states,  to  demonstrate  the  trends  observed  in  both  critical 
diameter  and  critical  initiating  pressure  with  increasing  charge 
homcgeneity. 

TNT  was  chosen  for  the  present  work  because  it  has  been  so 
widely  studied  that  more  data  are  available  for  it  than  for  any 
other  pure  explosive.  Moreover,  it  Is  a  castable  material  and, 

,  hence,  can  provide  an  almost  continuous  range  of  physical  hetero¬ 

geneity  between  the  highly  compacted  porous  charge  and  the  perfect 
single  crystal. 

THRESHOLD  FOR  SHOCK-TO- DETONATION  TRANSITION 

The  critical  initiating  pressure  Is  defined  as  the  minimum 
pressure  (P^)  required  to  initiate  detonation  of  the  explosive  in 
50#  of  the  trials.  It  defines  a  threshold  for  the  initiation  of 
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detonation  under  transient  conditions.  The  value  of  measured  in 
any  given  test  i3  determined  not  only  by  the  chemical  composition 
of  the  test  material  and  its  physical  state  but  also,  in  large  part, 
by  the  physical  dimensions  of  the  test.  These  dimensions  control 
the  times  at  which  rear  and  lateral  rarefactions  can  arrive  at  the 
shock  front  and  relieve  the  pressure.  Thus,  each  test  devised  to 
measure  is  apt  to  use  a  different  pressure  pulse  (pressure- time 
profile)  to  excite  reaction  of  the  explosive.  Different  profiles 
can  be  obtained  by  varying  the  composition  and  shape  of  donor 
explosives  in  the  gap,  booster,  and  wedge  tests  and  by  varying 
projectile  material,  shape,  and  velocity  in  the  projectile  Impact 
ana  wedge  experiments.  They  are  aloo  obtained  by  varying  the  diameter 
of  the  acceptor  explosive.  The  explosive  response  differs  as  the 
stimulus  differs,  and,  hence,  gives  different  P^  values  in  different 
tests.  In  other  words,  not  only  P^^  but  also  P^  vs  time  (t)  3hould 
be  specified  in  each  test. 

Eleven  years  ago  it  was  shown  theoretically  (l)  that  a  limiting 
pressure-time  (P-t)  curve  exists  for  a  50#  chance  of  build-up  to 
detonation  in  homogeneous  explosives;  the  limiting  curve  divides  the 
impulse  plane  into  a  detonation  region  and  a  failure  region.  Brown 
and  Whitbread  (2)  demonstrated  by  Impact  of  cylindrical  projectiles 
on  two  physically  heterogeneous  explosives  that  the  P^  determined 
must  have  a  minimum  duration  for  detonation  to  occur.  They  also 
showed  qualitatively  that,  at  P  >  P^  and  a  duration  shorter  than 
that  required  at  the  critical  pressure  level,  detonation  could  also 
be  initiated.  Since  then  there  has  teen  rather  general  agreement 
that  Initiation  of  detonation  must  be  the  result  of  the  pressure-time 
history  of  the  initiating  shock,  and  that  a  critical  limit  curve  in 
the  P-t  plane  must  exist  for  each  explosive. 

In  many  cases  of  explosive  loading  (e.g.,  a  gap  test),  maximum 
shock  pressure  can  be  determined  much  more  readily  than  the  pressure- 
time  profile.  Hence,  many  of  the  shock  sensitivity  measurements 
have  been  restricted  to  pressure  measurements  only.  However,  Walker 
and  Wasley  (3)  have  recently  combined  some  plate  projectile  impact 
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data  (for  which  pulse  duration  can  be  computed  from  plate  thickness) 
with  gap  test  measurements  about  which  some  information  on  the 
pressure-time  profile  was  available.  These  data  appear  to  define 
a  critical  ignition  energy  from  which  a  critical  pressure-time  limit 
curve  for  ignition  can  be  derived.  Walker  and  Wasley  also  discuss 
the  case  of  the  "sustained"  pulse,  i.e.,  one  of  duration  far  longer 
than  that  required  at  a  given  pressure  level  to  effect  ignition. 

They  believe  that  ignition  occurs  as  soon  as  the  critical  amount 
of  energy  has  been  transferred  to  the  test  explosive  from  the 
sustained  pulse. 

Some  of  the  gap  test  data  used  in  the  above  work  (3)  were 
obtained  with  the  calibrated  NOL  large  scale  gap  test  (LSGT)  (4). 
Because  it  will  be  used  for  data  in  this  work,  its  important 
characteristics  will  be  briefly  reviewed.  It  is  a  conventional 
gap  test  with  a  5  cm  diam  x  5  cm  long,  1.51  g/cc  tetryl  donor  and 
an  attenuator  gap  of  polymethyl  methacrylate  (PMMA)  or  its  equivalent. 
The  acceptor  explosive  can  be  tested  unconfined  or  in  the  moderate 
confinement  of  a  steel  sleeve  {3*65  cm  ID,  4.76  cm  OD).  A  mild 
steel  witness  plate  Is  used  in  both  cases,  and  the  gap  length  is 
varied  until  the  50#  value  is  found,  that  is,  the  attenuation  at 
which  a  hole  Is  punched  in  the  witness  plate  in  50#  of  the  trials. 

The  system  donor/gap  has  been  calibrated  (5)  to  give  shock  pressure 
as  a  function  of  gap  length.  This  50#  pressure  (P  )  at  the  end  of 
the  gap  can  be  converted  to  the  initiating  pressure  (P1)  transmitted 
to  the  explosive  by  use  of  the  Hugoniot  of  the  gap  material  (5) 
and  that  of  the  explosive,  e.g,,  cast  TNT  (6). 

The  standard  donor  of  the  LSGT  Is  approximately  point  initiated; 
hence,  the  detonation  front  in  the  donor  is  spherical  and  the  hydro- 
dynamic  flow  in  the  gap  and  in  the  acceptor  is  divergent.  The 
transmitted  pressure  falls  off  rapidly.  Prom  optical  observations 
of  the  tetryl/PMMA  system,  the  shock  pressure  decreases  to  half  Its 
initial  amplitude  in  about  one  psec  (7).  Walker  and  Wasley  (3) 
estimated  that  such  a  pulse  was  equivalent  to  a  square  pulse  of 
an  amplitude  of  0,9  P^  and  1.6  psec  duration.  Although  a  one¬ 
dimensional  hydrodynamic  flow  computation  on  the  tetryl/PMMA  system 
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produced  a  half-width  of  4  -jaec  (7),  analogous  two-dimensional 
computations  (8)  indicate  values  of  1.1  to  1.7  usee,  in  fair  agree¬ 
ment  with  the  original  estimate  from  experimental  work. 

Like  all  such  tests,  the  LSOT  exhibits  a  diameter  effect. 
Inasmuch  as  doubling  the  length  of  the  standard  booster  does  not 
affect  the  measured  pressure  (9)>  rear  rarefactions  play  no  part 
in  the  attenuation.  Thus,  the  observed,  as  well  as  the  computed, 
pulse  shape  can  be  attributed  to  the  effect  of  lateral  rarefactions 
on  the  initial  pulse.  Because  confinement  of  the  charge  delays  the 
arrival  of  lateral  rarefactions  at  the  charge  axis,  we  find,  for  the 
same  measured  P^,  that  the  diameter  of  the  unconflned  charge  i3 
approximately  twice  the  core  diameter  of  the  confined  charge,  i.e., 

76  mm  (See  Appendix).  This  value  of  the  equivalant  diameter  (dg) 
seems  applicable  to  a  number  of  different  explosives,  but  the  size 
of  the  change  in  measured  P^  with  confinement  depends  on  the  shock 
sensitivity  of  the  explosive.  The  change  is  greatest  for  the  least 
sensitive  materials  (See  Appendix). 

Both  pressed  and  cast  TNT  exhibit  the  usual  sensitivity 
phenomena  of  shocked  solid  heterogeneous  explosives  (10,11 ):  a 
critical  initiating  pressure  (P^  for  detonation,  breakout  of 
detonation  downstream  from  the  shocked  boundary,  and  a  consistent 
decrease  of  both  the  run  distance  and  delay  time  to  3teady-state 
detonation  with  increased  amplitude  of  the  applied  shock.  The 
run  distance  is  the  distance  from  the  plane  of  shock  entry  into 
the  explosive  to  the  plane  in  which  staady-state  detonation  is  first 
established.  The  delay  time  is  the  total  time  from  the  moment  of 
shock  entry  to  the  time  at  which  steady-state  detonation  begins. 

Liquid  TNT  and  presumably  single  crystal  TNT  behave  as  shocked 
homogeneous  explosives  (12).  They  too  exhibit  a  critical  initiation 
pressure,  but  detonation  occurs,  after  a  delay  (Induction)  time, 
at  the  plane  of  shock  entry.  The  detonation  wave  travels  through 
shocked  explosive  to  overtake  the  shock  wave.  There  is  subsequently 
a  short  period  of  overdrive  of  the  unshocked  explosive.  The  measured 
induction  time  decreases  with  increasing  pressure  of  the  Initiating 
shock.  It  has  been  shown  that  a  single  crystal  of  RDX  can  be 
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initiated  to  detonation  with  lower  strength  shocks  than  those 
required  for  homogeneous  initiation  (13)*  liquids  can  also  be 
initiated  at  similarly  lower  pressures,  e.g.  nitromethane  (14)  and 
nitroglycerine  (15).  These  appear  to  be  special  cases  in  which 
the  experimental  design  is  such  that  the  materials*  homogeneity 
can  be  destroyed  before  initiation  or  interaction  of  the  Initial 
shock  with  the  boundary  material  is  required  for  initiation.  Such 
special  cases  will  not  be  considered  further. 

LSCfT  shock  sensitivity  values  for  pressed  TNT  have  been  reported 
in  previous  work  (16)  as  Pg  vs  percent  theoretical  maximum  density 
{#  15®).  The  trend  was  the  usual  one  of  increasing  Pg  with  increas¬ 
ing  compaction.  The  average  particle  size  of  the  TNT  used  for  that 
work  was  about  150  u  with  a  maximum  particle  size  of  600  y.  To 
determine  the  effect  of  the  initial  particle  size  on  the  measured 
Pg,  a  special  batch  of  TNT  was  prepared  by  air  cooling  a  sprayed 
melt.  It  too  had  a  weight  mean  particle  size  of  150  p  by  sieve 
analysis.  It  was  separated  into  two  fractions:  that  passing  and 
that  retained  on  a  100  mesh  screen  (149  u  opening).  When  these 
materials  of  mean  particle  sizes  of  100  and  200  w,  respectively, 
were  compacted  isostatlcally  to  80#  15©  and  tested  In  the  standard 
confinement  of  the  LSGT  the  Pg  values  were  13.0  and  12.6  kbar 
respectively  as  compared  to  12.8  kbar  from  the  smoothed  curve  of  the 
earlier  work.  Although  these  differences  are  in  the  direction 
expected  from  test  results  on  other  explosives,  they  are  too  small 
to  be  significant  and  we  conclude  that  there  is  no  particle  size 
effect  evident  in  this  range  for  Pg  measured  on  pressed  TNT  In  the 
LSGT. 

We  have  measured  a  particle  size  effect  on  LSGT  values  for 
the  much  less  shock  sensitive  explosive,  nitroguanldine  (17).  In 
that  case  the  effect  was  small;  it  was  larger  at  lower  #  TMD  and 
disappeared  at  high  compactions.  The  coarser  material  showed  a 
slightly  greater  shook  sensitivity.  (TNT  was  tested  at  80#  TMD 
to  take  advantage  of  the  greater  effect  of  particle  size  on  Pg  at 
lower  compactions.)  In  gap  tests  of  smaller  effective  diameter 
than  the  LSGT,  the  trends  found  for  an  insensitive  explosive 
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(nitroguanldine)  have  be*n  confirmed  with  more  sensitive  materials, 
e.g,,  tetryl  and  PETN  (18). 

The  LSGT  measurements  on  TNT  have  been  collected  in  Table  1. 

Pg  values  have  been  converted  to  values  by  using  the  Hugoniots 
of  ?&MA  (5)  and  of  cast  TNT  (6)  except  for  the  one  low  density 
charge.  In  that  case  an  interpolation  was  made  between  the  Hugonlot 
for  TNT  (p0  *  1  g/cc)  (19)  and  that  at  higher  density  (6).  P^ 
values  at  d  >  3  d  (the  critical  diameter)  and,  hence  approximating 
the  infinite  diameter  value  are  given  for  each  charge.  A  single 
TNT  crystal  has  not  been  tested,  but  liquid  TNT  has  been  (12)  and 
3ur>plies  the  limiting  value  of  125  kbar.  [  It  was  assumed  that  the 
test  was  comparable  to  that  on  nitromethane  (12),  i.e.,  run  at  a 
diameter  of  about  8  inches.  Since  d  of  TNT  ({.)  at  8l  -  85°C  is 
reported  as  62  nun  (20)  and  68  mm  (21),  (d  /drt)  ~  3.  ]  Hie  initiating 
pressure  measured  lor  TNT  (-?)  is  a  limiting  value  for  the  crystal 
because,  from  the  Hugoniots  of  the  solid  (6)  and  of  the  liquid  TNT 
(22),  there  is  an  energy  Jump  across  a  125  kbar  shock  of  250  and 
360  cal/g,  respectively.  If  the  same  fraction  of  this  energy  is 
converted  to  heat  in  both  cases,  the  temperature  of  the  shocked 
solid  will  be  much  lower  than  that  of  the  shocked  liquid  because 
cv  (solid)  >  cy  (liquid)  (22).  Consequently  if  125  kbar  is  required 
for  homogeneous  initiation  of  the  liquid,  a  higher  pressure  would 
be  required  to  initiate  the  solid.  The  validity  of  this  conclusion 
is  also  supported  by  a  recent  calculation  of  lor  a  single  TNT 
crystal.  Using  a  measured  adiabatic  explosion  time,  Arrhenius 
kinetics,  and  an  assumed  equation  of  state  for  the  solid, 
Voskoboynikov  et  al  (23)  computed  a  critical  initiating  pressure 
of  200  kbar  for  the  perfectly  homogeneous  solid  TNT. 

"Values  for  the  two  cold  pressed  charges  of  Table  1  are  smoothed 
values  from  previous  work  (16).  The  values  for  the  confined  charges 
of  two  of  the  castings  have  been  deduced  from  measurements  on  bare 
charges,  as  described  in  the  Appendix.  All  values  have  been  arranged 
in  order  of  increasing  physical  homogeneity  of  the  charge.  Judged 
from  a  knowledge  of  its  method  of  preparation.  Thus,  a  rapid  cooling 
and  solidification  of  TNT  results  in  a  casting  containing  many  small 
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crystals;  it  should  be  the  most  heterogeneous  of  the  castings  although 
more  homogeneous  than  any  pressed  charge.  A  vacuum  casting  should 
be  more  homogeneous  than  a  creamed  casting  both  because  of  the  evacu¬ 
ation  of  any  trapped  air  and  because  cooling  under  vacuum  would  be 
somewhat  slower  than  the  air  cooling  of  the  creamed  casting.  Casting 
Wo.  4  was  obtained  with  a  heated  mold  and  steam  fingers  which  supplied 
some  heat  to  the  interior  of  the  casting  throughout  the  solidification 
process.  These  charges  contained  three  thin  metal  wires  (one  along 
the  axis,  the  others  parallel  to  it  with  a  separation  of  3  mr.)  which 
served  to  conduct  the  heat.  The  charges  so  prepared  contained  a 
number  of  large  TNT  crystals  and  most  closely  approximated  the  case  of 
the  single  crystal. 

Although  a  single  operator  might  learn  to  reproduce  any  of  these 

castings  fairly  well,  the  variation  from  operator  to  operator  and  from 

day  to  day  can  be  large  despite  written  schedules  and  procedures.  For 

example,  eight  different  preparations  of  creamed  cast  TNT  showed  a 

variation  of  26.7  to  46.4  kbar  in  ?g.  If  the  two  extreme  values  are 

discarded,  six  charges  showed  Pg  values  of  33  to  40  kbar.  In  view  of 

the  sensitivity  of  P_  to  small  changes  in  mold  temperature  and  cooling 

© 

rate,  this  range  is  probably  as  small  as  could  be  expected.  The  two 
values  measured  for  Casting  2  of  Table  1  (confined  and  unconfined) 
were  obtained  from  charges  made  from  the  same  lot  of  TOT.*,  and  the  P 
measured  under  confinement  is  the  same  as  the  mean  value  obtained  for 
either  the  group  of  eight  or  that  of  six  creamed  cast  charges. 

The  data  of  Table  1  show  two  things  very  clearly:  (1)  there  is 
a  continuous  variation  in  values  from  the  most  to  the  least  hetero¬ 
geneous  charge,  and  (2)  the  required  initiating  pressure  increases 
monotonically  with  increasing  physical  homogeneity  of  the  charge. 
Although  a  smooth  trend  in  the  range  of  cold  pressed  charges  is 
commonly  found,  earlier  work  (16)  suggested  a  discontinuity  between 
the  initiating  pressure  required  for  the  highest  density  pressed 
charge  and  that  for  a  single  crystal  (P.DX  or  PETN) .  If  the  values  for 
the  cast  charges  were  omitted  from  Table  1,  a  similar  apparent  Jump 
in  would  occur  for  TNT.  A  castable  material  can  be  handled  to 
produce  charges  of  homogeneity  bridging  the  difference  between  that  of 
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the  highest  density  pressed  charge  and  of  the  perfect  crystal.  A 
comparable  treatment  for  non-castable  explosives  such  as  RDX  would  be 
the  use  of  a  solvent,  i.e.,  solution  in  place  of  melting.  In  both 
cases,  a  continuous  variation  in  homogeneity  and  hence  in  P.^  can 
probably  be  obtained. 

Finally  it  should  be  noted  that  pQ  or  %  TMD  which,  together  with 
the  particle  size  of  the  explosive,  seems  an  important  variable  in 
determining  for  pressed  charges,  has  little  effect  on  P^  in  the 
range  of  the  castings.  The  controlling  variable  for  the  entire  range, 
as  mentioned  above,  is  the  degree  of  physical  homogeneity  of  the  charge. 
Unfortunately,  there  is  as  yet  no  quantitative  measure  of  this  charge 
characteristic;  however,  for  a  chemically  homogeneous  compound  such  as 
TNT,  degree  of  homogeneity  would  certainly  encompass  the  bulk  density 
(oQ)  with  an  appropriate  weighting  for  its  greater  apparent  Influence 
on  the  behavior  of  the  more  heterogeneous  charges.  Qualitatively  it 
Is  clear  that  the  shock  sensitivity,  as  measured  by  P decreases 
montonically  with  increasing  physical  homogeneity. 

THRESHOLD  FOR  PROPAGATION  OF  DETONATION 

The  failure  or  critical  diameter  (d  )  determines  the  charge  size 
at  which  there  exists  a  threshold  for  the  propagation  of  detonation 
under  steady  state  conditions.  By  definition,  dc  is  that  diameter  of 
a  cylindrical  charge  at  or  above  which  detonation  propagates  and  below 
which  it  fails.  The  existence  of  a  failure  diameter  is  the  result  of 
two  dimensional  effects.  In  accord  with  this  and  with  the  shock-to- 
deton&tion  transitions  discussed  in  the  previous  section,  dc  defines 
that  charge  size  at  which  the  steady-state  detonation  wave  has  been 
attenuated  from  its  infinite  diameter  value  by  lateral  rarefactions 
until  it  is  Just  critical  for  initiation.  In  other  words,  at  dQ  the 
pressure  pulse  between  the  von  Neumann  shock  front  and  the  C-J  plane 
is  Just  critical  for  initiating  detonation  after  a  run  length  equal  to 
the  reaction  zone  length  and  a  total  delay  time  equal  to  the  reaction 
time.  Thus  at  d  ,  the  transient  initiation  phenomena  just  fit  into 
the  reaction  zone  at  the  threshold  of  steady-state  propagation.  This 
is  the  rather  complex  relationship  between  the  explosive  behavior  at 
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these  different  thresholds,  each  of  which  is  a  different  facet  of 
shock  sensitivity.  A  quantitative  illustration  of  this  proposed 
relationship  will  be  given  later. 

Available  Data  on  d„  of  TNT 

c _ 

In  the  range  of  cold  pressed  charges,  the  variation  of  the  criti¬ 
cal  diameter  with  the  charge  porosity  has  served  as  the  basis  of  a 
convenient  classification  of  explosives  into  two  groups  (16).  Group  1, 
typified  by  TNT,  exhibits  decreasing  dQ  with  increasing  compaction 
(increasing  %  TMD) .  Group  2,  typified  by  ammonium  perchlorate  (AP) , 
exhibits  the  reverse  trend.  As  a  result  of  the  different  failure 
behavior,  the  two  groups  also  exhibit  differences  in  their  detonation 
velocity  (D)  vs  pQ  curves  at  finite  diameters.  Members  of  Group  1 
have  linear  D  vs  pQ  curves  whereas  those  of  Group  2  can  exhibit  D  vs 
pQ  curves  with  a  maximum  in  D.  In  later  work  (17)  we  found  that  the 
division  between  the  two  groups  was  not  as  sharp  as  we  first  thought, 
that  some  materials,  e.g.,  nitroguanidine,  have  a  U  shaped  dc  vs  pQ 
curve  and  exhibit  Group  1  patterns  at  lower  %  TMD,  Group  2  patterns 
at  higher  %  TMD.  It  is  probable  that  any  pressed  explosive  can  behave 
in  this  way  provided  the  necessary  range  in  degree  of  compaction  can 
be  achieved  experimentally. 

Most  of  the  d  measurements  on  pressed  TNT  (porous,  granular 
compacts)  have  been  made  by  Russian  investigators.  Pig.  1  shows  some 
of  these  data  plotted  dc  vs  pQ.  The  two  solid  lines  are  values 
reported  by  Bobolev  (24)  for  fine  and  coarse  TNT.  These  curves  show 
the  Group  1  behavior  described  above  and  also  demonstrate  the  effect 
of  particle  size  on  critical  diameter.  In  the  oQ  range  of  O.85  to 
1*23  g/cc,  d  is  decreased  by  a  factor  of  2  to  3  for  a  decrease  in  the 

V 

Initial  average  particle  diameter  (increase  in  surface  area)  of  about 
a  factor  of  4.  It  should  be  noted,  however,  that  the  curve  for  the 
fine  TNT  ends  at  75$,  probably  because  of  pressing  difficulties.  If 
experimental  data  could  be  obtained  to  extend  that  curve,  it  should 
approach  the  curve  for  the  initially  coarser  TNT  even  as  the  two 
vs  pc  curves  become  coincident  at  high  compactions.  That  this  is 
probably  the  case  is  indicated  by  a  recent  determination  (25)  on  fine 
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CRITICAL  DIAMETER  dr  (mm) 


FIG.  1  FAILURE  LIMIT  CURVES  FOR  TNT  IN  DIAMETER  -  DENSITY  PLANE 

( - REF  24;  x  REF  25;  •  DETONATION  O  FA  I  LURE,  REF  28; 

- CURVE  CHOSEN  FOR  REF  28  DATA;  -  —  --INTERPOLATION) 
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(20  -  70u)  TNT  at  97#  TVS).  In  this  work  a  critical  slab  thickness 
(h)  was  measured;  it  and  the  approximate  d  ~  2h  are  plotted  as 

C  — 

crosses  in  Figure  1.  Either  value  indicates  that  the  fine  TNT 
critical  curve  will  approach  or  go  through  a  minimum  and  that  at 
high  compaction  there  will  be  no  significant  difference  in  d  of 

V 

the  initially  fine  and  initially  coarse  TNT. 

Although  there  is  only  one  set  of  data  for  the  fine  TNT,  there 
are  several  for  the  coarser  material.  Dremin  et  al  (26)  quote  limit 
data  from  Apin  and  Stesik  (27)  for  100  u  TNT.  That  failure  limit 
curve,  not  shown  in  Figure  1,  lies  0.8  to  0.2mm  below  Bobolev's 
curve  for  70  -  200u  TNT,  about  where  it  should  be  for  a  small 
decrease  in  the  average  particle  size.  The  most  recent  data  are 
those  of  Stesik  and  Akimova  (28);  their  limiting  shots  for  failure 
and  detonation  are  plotted  in  Figure  1.  Over  the  range  1.0  to  1.62 
g/cc  they  agree  quite  well  with  the  data  of  reference  24;  at  0.8  g/cc 
there  is  a  large  disagreement.  The  close  agreement  at  higher  densi¬ 
ties  makes  the  size  of  the  discrepancy  between  the  most  porous 
charges  seem  Improbable.  Examination  of  the  D  vs  reciprocal  diameter 
(d-1)  curve  (28)  for  the  charges  at  pe  =  0.8  g/cc  shows  a  sharp  change 
in  slope  at  d  -  10mm.  For  this  reason,  d  (0.8  g/cc)  =  10mm  is  used 
here  instead  of  the  lower  value  selected  in  reference  28.  (At  d  < 
dc,  it  is  often  difficult  to  distinguish  an  overdriven  shock  from  a 
steady-state  detonation  in  a  low  density  charge  of  moderate  length.) 
With  this  change,  the  reference  28  data  are  in  far  better  agreement 
with  data  from  both  reference  27  and  reference  24,  as  shown  in 
Figure  1. 

The  possible  forms  of  the  lower  curve  of  Figure  1  at  higher 
density  strongly  suggests  that  TNT  can  be  put  in  a  physical  state  in 
which  it  will  exhibit  a  U  shaped  d  vs  p0  curve .  This  possibility  has 
been  recently  confirmed  with  measurements  on  hot  pressed  (72°  -  76°C) 
TNT.  In  the  narrow  range  of  1.60  to  1.65  g/cc  in  p0  and  8  to  13mm 
in  d,  these  charges  exhibit  Group  2  patterns  for  both  D  vs  p0  and 
dc  vs  o0  curves  (29)  i.e.,  a  maximum  D  in  the  D  vs  p0  curve  and  dQ 
increasing  with  increasing  p0 .  The  data  are  plotted  in  Figure  2. 
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The  above  measurements  on  pressed  charges  are  summarized  in 
Table  2  which  also  includes  available  data  for  cast  TNT.  The  data 
have  been  arranged  In  order  of  increasing  physical  homogeneity  of  the 
charge,  as  were  those  of  Table  1.  As  in  the  case  of  P^,  there  is  a 
continuous  variation  in  d  from  the  most  to  the  least  heterogeneous 
charge.  The  gap  between  the  values  for  the  highest  density  pressed 
charge  and  the  single  crystal  (with  d  greater  than  that  of  Casting  4) 
is  smoothly  bridged  by  charges  produced  by  the  different  methods  of 
preparation.  Similarly,  the  gap  between  the  cold  pressings  (hetero¬ 
geneous  and  permeable  to  fluids)  and  the  castings  (approximately 
homogeneous  and  Impermeable)  can  be  bridged  by  hot  pressing.  Pure 
TNT  melts  at  8l°C;  hence,  pressing  at  72  -  76°C  will  cause  some  local 
melting  and  flow  so  that  the  charges  contain  both  permeable  and 
impermeable  regions.  In  the  region  of  cold  pressed  charges,  d 

c 

decreases  with  increasing  degree  of  homogeneity  (a  Group  1  pattern), 

but  this  trend  is  reversed  in  the  region  of  hot  pressing.  Thereafter, 

dc,  like  P^,  increases  with  increasing  homogeneity.  The  overall 

trend  strongly  suggests  that  permeability,  permitting  convective  heat 

transfer  by  flow  of  hot  gas  products,  has  a  role  in  the  initiation  and 

propagation  of  detonation.  Finally  the  data  of  Table  2  show  that 

the  effect  of  the  initial  particle  size  on  d  is  appreciable  at  low 

c 

%  TMD  and  negligible  at  high  %  IMP,  as  it  was  for 
Other  Information  about  Critical  Conditions 


More  information  can  be  obtained  about  the  failure  conditions 
if  other  measurements  are  made  in  addition  to  dc.  In  particular, 
from  the  detonation  velocity  at  d_  (D  ),  an  estimate  of  the  detonation 

C  C 

pressure  at  critical  conditions  Is 


pj=  s  p°  V/4  U) 

A  measure  of  the  reaction  zone  length  (x)  leads  to  an  estimate  of 
reaction  time  (T).  Since  x  =  (D  -  u)t  where  u  is  the  average  particle 
velocity  in  the  reaction  zone,  the  estimate  "u  =  1.2  Uj  leads  to  x^^  = 
f(y)  where  f(?)  is  a  slowly  varying  function  of  the  adiabatic 

exponent  gamma;  f(y)  can  be  approximated,  as  in  reference  (33), 
by  2/3  to  give 
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Tj,  «  1.5  (2) 

Equation  (2)  gives  the  approximate  value  for  the  Ideal  or  Infinite 
diameter  detonation  reaction  time  because  measurements  or  estimates 
of  x  are  made  at  conditions  approaching  ideal.  For  near  voidless 
materials,  however,  Dc  =  and  should  be  a  good  approximation 
for  t  .  Finally  an  estimate  of  the  energy  required  for  propagation 

C 

at  dfi  can  be  obtained  from 

E  (cal/cm3 )  -  P® t/p0  U  (3) 

Where  U  Is  the  shock  velocity  of  amplitude  P  in  the  unreacting 
explosive.  Walker  and  Wasley  (3)  derived  Equation  (3)  for  a  square 
pressure  pulse  of  amplitude  P  and  duration  t  created  by  high  velocity 
impact  of  an  explosive  projectile  against  an  explosive  target. 

If  P*  T»  Po  >  and  U  are  in  units  of  kbar,  psec,  g/cm3,  and 
mm/^sec,  respectively,  the  conversion  factor  is  0.239.  With  the 
hydrodynamic  relation,  P  =  p0Uu,  Equation  (3)  can  be  converted  to 

E  =  PTu  (4) 

Where  u  is  particle  veloc5ty.  For  the  units  of  kbar,  usee  and 
mm/usec,  a  factor  of  2.39  gives  E  in  cal/cm3.  To  obtain  the  energy 
value  at  d^,  where  U  =  Dc,  we  approximate  P  ~  Pjc  and  tc  =  t^.  The 
result  is  a  low  estimate  because  the  square  pressure  pulse  equivalent 
to  the  pulse  between  the  von  Neumann  and  C-J  planes  should  have  an 
amplitude  greater  than  PJc  (Tj,  is  probably  less  than  rc  as  well). 

Steslk  and  Akimova  (28)  reported,  in  addition  to  the  d  values 

C 

already  used,  a  number  of  D  vs  d  measurements  and  also  derived 
reaction  zone  length  data.  Their  data  and  results  computed  from  them 
with  Equations  (l)  to  (4)  are  summarized  in  Table  3.  All  the 
critical  parameters  of  the  pressed  charges  show  a  smooth  variation 
with  p0  or  %  TMD,  but  the  only  linear  curves  are  dQ  vs  $  TMD  and 
log  Pjc  vs  $  TMD.  The  rest  are,  like  (D(J/D1)  vs  $  TMD,  concave 
upward.  The  tabulated  values  of  (^/D^)  show  the  low  values  (ca.  0.6) 
typical  of  loose  powder  -barges  and  the  high  values  (ca.  0.95)  typical  of 
cast  or  pressed  charges  at  95$  or  more  of  their  voidless  density 
(see  Figure  3)»  The  detonation  pressure  at  the  failure  limit  is. 
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of  course,  a  smaller  fraction  of  the  corresponding  infinite  diameter 
value  than  is  DQ  because,  by  the  approximation  of  Equation  (1), 

V»n  -  <W- 

ISiere  is  at  present  some  confusion  about  measuring  detonation 
reaction  zone  length  and  reaction  time.  Table  3  gives  reaction 
zone  lengths  (x^)  and  reaction  time  (t^)  computed  from  Eyring's 
curved  front  theory.  Stesik  and  Akimova  (28)  obtained  an  experi¬ 
mental  value  of  the  zone  length  on  the  basis  of  calibrating  their 
system  with  one  series  of  free  surface  velocity  measurements  on 
explosively  driven  plates  of  different  thicknesses.  They  also 
computed  from  these  results  a  detonation  reaction  time  (t^)  which 
is  presented  for  comparison  in  Table  3.  The  times  t  were  chosen 
for  the  present  treatment  because  they  produce  a  linear  curve, 
log  t  vs  %  TMD,  which  is  also  a  fairly  good  approximation  to  the 
reference  28  data.  Both  t  and  t  are  of  the  same  order  of  magnitude 
as  the  more  recently  reported  direct  measurements  of  reaction  time 
(36).  The  reference  36  data  lie  between  the  t  and  t  values  of 
Table  3  at  pa  =0.8  to  1.0  g/cc  and  above  them  at  1.4  g/cc,  as 
Figure  4  shows.  No  comparison  is  made  for  results  on  cast  TNT 
since  there  is  no  way  to  be  sure  that  the  castings  of  references 
28  and  36  are  comparable. 

The  last  column  of  Table  3  contains  the  values  of  energy 
per  unit  area  (E)  delivered  to  the  explosive  to  sustain  propagation 
at  the  critical  conditions.  In  view  of  the  approximations  used, 

E  might  be  constant  over  the  range  of  the  pressed  charges  but  it 
is  certainly  lower  than  the  value  required  for  a  cast  charge.  It 
should  be  noted  that  if  Equation  (3)  is  converted  to  energy  per 
gram  instead  of  energy  per  unit  area,  it  is  equal  to  Ujc3  and  shows 
a  smooth  variation  with  %  1MD  as  do  the  otl  er  detonation  parameters 
at  d  »  dc. 

Finally  the  data  of  Table  3  can  be  used  to  construct  a  limit 
curve  in  the  pressure-time  plane  for  pressed  TNT  charges.  This  is 
shown  in  Figure  5  where  the  curve  divides  the  plane  into  a  super- 
and  a  subcritical  region. 
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FIG.  4  COMPARISON  OF  TNT  REACTION  TIMES  [o  r,  TABLE  3  ;  x  t,  TABLE  3  ;  A  REF  36] 
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DISCUSSION 

The  data  gathered  for  TNT  have  ahown  that  both  the  critical 
Initiating  pressure  F,  and  the  critical  diameter  d.  are  determined 

X  v 

in  large  part  by  the  degree  of  physical  homogeneity  of  the  charge. 

It  is  to  be  expected,  therefore,  that  these  two  threshold  values 
will  be  related  to  each  other,  and  the  plot  of  dc  vs  shown  in 
Figure  6  affirms  such  a  relationship  by  showing  a  smooth  carve 
over  the  entire  experimental  range.  Oliis  curve  also  serves  as  a 
graphic  summary  of  much  of  our  present  information  on  the  critical 
behavior  of  an  organic  explosive. 

Although  there  is  no  quantitative  measure  of  the  degree  of 
physical  homogeneity,  we  know  that  it  varies  from  very  low  at  the 
extreme  left  of  Figure  6  (high  porosity,  cold  pressed  charges)  to 
very  high  at  the  extreme  right  (100#  when  the  extrapolated  curve 
reaches  the  perfect  single  crystal).  Ofcese  extremes  are  indicated 
on  the  figure.  Between  these  extremes,  shows  a  mono tonic  Increase 
with  increasing  homogeneity  and  no  particle  size  effect  in  the  range 
of  cold  pressed  charges.  (That  is,  none  can  be  detected  down  to  80# 
THU  in  the  LSOT,  Other  tests  do  show  a  small  effect  for  other 
explosives.)  At  the  extreme  right  (single  crystal)  ignition  and 
initiation  must  be  by  a  homogeneous  mechanism  since  no  other  is 
available.  At  the  extreme  left,  it  is  generally  agreed  that  initia¬ 
tion  must  be  by  a  hot  spot  mechanism,  i.e.,  surface  reaction  at 
areas  of  energy  concentration.  This  heterogeneous  mechanism  is 
accepted  because  the  measured  initiating  pressure  is  far  too  low 
to  effect  by  shock  compression  a  significant  temperature  rise 
throughout  the  bulk  of  the  material.  As  the  charges  are  changed 
from  heterogeneous  to  homogeneous,  the  ignition  mechanism  also 
changes;  it  follows  that,  in  some  intermediate  region,  initiat^n 
must  be  by  both  heterogeneous  and  homogeneous  mechanisms  acting 
simultaneously.  The  effect  of  degree  of  homogeneity  on  is 
attributed  to  its  reflection  of  the  concentration,  distribution, 
and  size  of  hot  spot  sites. 
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VARIATION  OF  CRITICAL  DIAMETER  WITH  CRITICAL  INITIATING  PRESSURE 

FOR  SEVERAL  FORMS  OF  TNT 

(COLD  PRESSED  CHARGES:  □  COARSE,  ■  FINE; 

HOT  PRESSED  CHARGES - ;  CAST  CHARGES  O  ;  EXTRAPOLATION - 
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In  the  region  of  cold  pressed  charges,  a  particle  size  effect  on 
dc  has  been  indicated.  The  Group  1  behavior  of  decreasing  dc  with 
increasing  %  TliD,  shown  by  the  coarse  TUT,  was  previously  attributed 
to  increasing  declination  of  the  heterogeneous  reaction  by  the  homo¬ 
geneous  (16).  Thus,  at  higher  density  and  higher  reaction  pressure, 
the  dominance  of  bulk  over  surface  reaction  could  result  in  shorter 
reaction  times  and  reaction  zone  lengths  which,  in  turn,  would  result 
in  lower  d  values.  This  still  seems  a  reasonable  suggestion  because 

V 

detonation  pressures  of  120-180  kbar  (TAble  3  and  Fig.  5)  are  quite 
high  enough  to  cause  appreciable  bulk  heating  during  shock  compression. 
But  the  apparent  reversal  in  trend  shewn  by  the  fine  TNT  and  the 
obvious  reversal  in  trend  (dc  vs  or  dc  vs  degree  of  homogeneity) 
in  the  hot  pressed  charge  range,  as  compared  to  the  cold -pressed, 
suggests  that  another  factor  should  be  considered.  Because  the 
behavior  of  the  hot  pressed  charges  bridges  that  of  the  cold-pressed 
and  cast  which  are,  respectively,  permeable  and  impermeable,  the  new 
factor  might  well  be  the  flow  of  detonation  products  through  the 
charge.  If  so,  convective  heat  transfer  from  such  flow  must  play  a 
role  such  that  initiation  is  easier  to  achieve  when  the  flow  is  optimum 
in  permeable  charges.  Charges  made  of  initially  coarse  explosive  are 
more  permeable  than  those  made  of  fine  (at  the  same  %  TMD)  until  the 
compaction  eliminates  any  difference.  Hence  it  may  well  be  a  difference 
in  permeability  rather  than  in  total  surface  area  that  is  responsible 
for  the  particle  size  effect  found  in  the  cold -pressed  Tiff  charges. 

This  suggestion  is  strengthened  by  the  apparent  importance  of 
permeability  in  the  transition  from  burning  to  detonation  in  granular 
compacts (37) .  Although  permeability  was  not  the  only  factor  affecting 
the  ease  of  transition  (for  example,  surface  to  volume  ratio  of  the 
individual  particles  and  energy  produced  by  chemical  reaction  were  also 
considered),  Griffiths  and  Groocock  (37)  f~und  it  a  very  important 
one.  They  proposed  that  the  high  velocity  combustion  necessary  for 
transition  was  possible  because  of  a  convective  transfer  of  energy 
from  hot  reaction  gases  to  solid  unreacted  particles.  Moreover,  their 
curve  of  the  length  of  burning  before  detonation  vs  permeability  of 
HMX  charges  pressed  to  different  densities  was  very  similar  to  the 
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curves  they  obtained  by  pressing  to  a  constant  density  and  varying  the 
permeability  by  varying  the  initial  particle  size.  A  shock-to- 
detonation  transition  is  a  very  rapid  burning -to -detonation  transition 
in  which  experimental  time  resolution  is  difficult.  Mechanisms 
operative  in  more  easily  resolved  transitions (37)  may  also  be  opera¬ 
tive  here. 

In  this  connection  it  should  be  noted  that  the  phenomenon  of 
"overshoot"  observed  in  the  shock  initiation  of  detonation  in  cold 
pressed  charges  (10,  38,  39)  is  not  observed  in  hot  pressed  charges 
(4o)  or  in  cast  charges  (10).  The  extent  of  overshoot  in  cold  pressed 
charges  is  greatest  for  the  most  porous  and  least  for  the  most  compact 
(10).  By  arguments  similar  to  those  already  used,  permeability  of  the 
charge  may  be  responsible  for  the  appearance  of  overshoot. 

The  details  of  exactly  hew  a  forced  convective  flow  contributes 
to  the  chock-to-detonation  transition  have  not  yet  been  determined. 

The  initial  flow  would  be  expected  to  be  slower  than  the  supersonic 
progress  of  the  initial  shock  front.  Hence  any  reaction  it  induced 
would  occur  after  the  passage  of  the  shock  although  subsequent  growth 
of  that  reaction  could  result  in  producing  a  front  which  could  over¬ 
take  the  shock.  This  sequence  of  events  could  be  responsible  for  the 
observed  decrease  of  pressure  in  the  acceptor  with  the  length  of  the 
shock  path  into  the  explosive.  The  initial  decrease  of  pressure  is 
followed  by  an  increase  (at  greater  path  lengths)  when  a  successful 
transition  to  detonation  occurs  (4l) . 

Figure  6  summarizes  the  large  effects  of  the  physical  homogeneity 
on  the  propagation  and  initiation  thresholds  of  a  pure,  chemically 
homogeneous  explosive.  Also  indicated  are  the  suggested  mechanisms 
associated  with  the  different  ranges.  The  complexity  of  Fig.  6  illus¬ 
trates  the  impossibility  of  creating  a  single  satisfactory  model  for 
a  theoretical  treatment  of  either  failure  or  shock-to-detonation 
transition  over  the  entire  range  of  variation  in  initial  physical 
state.  At  the  same  time,  however.  Fig.  6  shows  that  there  are  no 
discontinuous  changes  in  the  threshold  values  when  the  preparation  can 
be  varied  to  give  charges  of  different  degrees  of  homogeneity  so  that 
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they  bridge  the  rather  large  gap  between  the  highest  density  pressed 
charge  and  the  perfect  single  crystal. 

Although  Fig.  6  demonstrates  the  continuum  of  threshold  values, 
it  does  not  provide  a  quantitative  example  of  the  relationship  between 
conditions  at  these  two  types  of  threshold  i.e.,  of  the  transient 
initiation  phenomena  occurring  within  the  steady  state  reaction  zone 
at  the  threshold  of  propagation.  We  can  do  this  in  two  ways:  (a)  plot 
the  two  sets  of  pres sure -time  data  (those  at  dc  and  those  fron  the 
LSGT)  on  a  log-log  scale,  as  suggested  by  Fig.  5,  and  approximate  the 
limit  curve  by  a  straight  line  or  (b)  use  the  concept  of  a  critical 
energy  for  initiation  of  detonation,  compute  that  energy  from  the 
measurements  at  dc  with  Eq.  (3)  ana  use  it  in  Eq.  (4)  to  define  the 
limit  curve.  In  either  case,  we  must  know  the  Hugoniot  (P  vs  u)  of 
the  non-reacting  explosive  or  approximate  it  in  some  way.  There  are 
several  Hugoniot s  for  high  density  TUT  which  can  be  used  for  the  near 
voidless  TOT,  pressed  or  cast,  but  there  is  only  one  for  a  very  porous 
(1.0  g/cc)  TOT  (19).  There  is  some  question  about  the  possibility  of 
shocking  such  a  porous  material  without  inducing  reaction,  i.e.,  the 
Hugoniot  data  may  be  those  of  a  partially  reacting  material  rather 
than  the  desired  data  for  the  non-reacting  material.  Moreover,  this 
very  porous  charge  is  in  the  region  where  the  effect  of  the  initial 
particle  size  is  greatest  (see  Fig.  6)  and  where  we  lack  the  necessary 
critical  data  on  charges  prepared  from  fine  TOT.  For  both  of  these 
reasons,  we  chose  to  use  pressed  TOT  of  pQ  =  1.62  g/cc. 

Table  3  shows  that  the  "critical  energy"  for  initiation  of  the 
high  density  compact  is  32  cal/cm  .  This  is  remarkably  and  probably 

p 

fortuitously  close  to  the  energy  limit  of  33  cal/cm  Walker  and  Wasley 
report  for  ignition  of  pressed  TOT  at  1.64  g/cc.  Our  value  is  low 
because  of  approximating  the  critical  pressure  pulse  with  a  square  one 
of  amplitude  Pjc  (a  minimum  rather  than  average  value)  but  it  does 
3eem  to  be  of  the  right  order  of  magnitude.  It  is  combined  with  the 
Hugoniot  data  (6)  in  Eq.  (4)  to  generate  the  limit  curve  of  Fig.  7. 

The  curve  is  slightly  concave  upward,  but  can  be  fairly  well  approxi¬ 
mated  by  a  straight  line.  Moreover  the  estimated  square  pulse 
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equivalent  of  the  intitlating  pulse  used  in  the  LSGT  (0.9  1.6  us®c) 

falls  on  the  limit  curve  as  shown,  well  within  errors  tc  be  expected 
from  those  introduced  experimentally  as  well  as  by  the  approximations 
used.  The  limit  curve  divides  the  P-t  plane  into  a  possible  Initiation 
and  a  failure  region  (above  and  below  the  curve,  respectively).  The 
curve  originates  at  the  P-t  conditions  found  at  the  critical  diameter 
threshold  and  runs  through  the  continuum  of  (P,t)  values  which  in- 
•eludes  all  the  threshold  values  for  initiation  measured  in  the  various 
gap,  wedge,  projectile,  and  booster  tests.  The  conditions  for  the 
large  scale  gap  test  fall  on  the  curve  toward  its  lower  end;  the  curve 
probably  terminates  shortly  beyond  that  point  since  there  is  probably 
an  upper  time  limit  in  which  the  critical  energy  must  be  delivered  in 
order  to  initiate  detonation.  This  particular  example  provides, 
therefore,  a  consistent  and  satisfying  illustration  of  the  conceptual 
relationship  between  conditions  at  the  threshold  for  propagation  and 
the  numerous  thresholds  for  initiation.  It  also  supports  the  concept 
of  a  critical  energy  of  Initiation,  but  further  tests  with  improved 
data  should  certainly  be  made. 

In  principle,  a  limit  curve  analogous  to  that  for  the  pressed 
TlfT  can  be  constructed  for  the  cast  TlfT.  In  fact,  energies  computed 
for  the  cast  charge  (rapidly  cooled)  of  Table  3  from  critical  diameter 
data  and  from  comparable  Pi  data  (Table  1)  differ  by  a  factor  of  about 
two.  This  discrepancy  is  attributed  to  error  in  the  reaction  time 
used  for  the  cast  TNT.  If  we  assume  that  the  highest  density  pressed 
charge  and  all  the  castings  have  the  same  values  Dc  and  Pjc,  as  the 
Table  3  data  suggest,  we  can  estimate  critical  energies  and  reaction 
times  from  Eq.  (3)  and  the  P^  data  of  Table  1.  With  the  LSGT  condi¬ 
tions  of  P  =  0.9  P^,  t  =*  1.6  usee,  Eq.  (3)  indicates  relative  reaction 
times  at  d  of  1,  1.6,  2.6,  and  5.9, respectively,  for  the  pressed, 

V 

rapidly  cooled  cast,  creamed  cast,  and  vacuimi  cast  charges.  This  is 
a  reasonable  trend  and  gives  an  increase  in  critical  energy  from  the 
high  density  compact  to  Casting  1  of  1.6  instead  of  4.1  times  as 
shown  in  Table  3. 
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SUMMARY 

The  various  aspects  of  shock  sensitivity  have  been  illustrated 
with  data  obtained  for  TNT  charges  of  varying  degrees  of  physical 
homogeneity.  Beth  the  critical  diameter  for  the  propagation  of  deto¬ 
nation  (d  )  and  the  critical  initiating  pressure  measured  with  the 

V 

LSGT  (P^)  show  continuous  change  with  charge  homogeneity.  increases 
(ease  of  initiation  decreases)  monotonically  with  Increasing  charge 
homogeneity.  The  same  trend  was  observed  for  d  in  the  region  of  cast 

w 

charges,  but  the  reverse  trend  occurs  in  the  region  of  cold  pressed 
charges.  The  reversal  occurs  in  che  region  of  hot  pressed  charges 
which  bridges  the  cold-pressed,  permeable  charges  and  the  cast, 
impermeable  ones.  The  reversal  is  associated  with  the  effect  of 
permeability  on  the  process  of  initiating  detonation. 

It  is  proposed  that  any  charge  will  exhibit  a  critical  curve  in 
the  pressure-time  plane  and  that  such  a  curve  will  start  at  a  point 
corresponding  to  the  detonation  pressure  and  reaction  time  at  d  .  The 

v 

limit  curve  will  then  run  through  a  continuum  of  (P,t)  values  corres¬ 
ponding  to  the  various  gap,  booster,  wedge,  and  projectile  test  values 
for  sensitivity.  In  an  example  for  pressed  TNT  at  p0  =  1.62  g/cc, 

P^  from  the  LSGT  fell  on  such  a  curve  computed  from  a  critical  energy 
of  initiation;  the  energy  was  evaluated  from  data  obtained  at  dc.  This 
example  supports  the  proposed  model  in  which  the  transient  Initiation 
phenomena  occur  within  the  reaction  zone  at  the  threshold  of  steady 
state  propagation.  It  also  supports  the  concept  of  critical  energy 
for  initiation  of  detonation.  Further  tests  of  both  concepts  should 
be  made  with  improved  data. 
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APPENDIX 

Equivalent  Diameter  for  ROL  LSQT  Confinement 

In  some  of  the  earliest  work  on  quantitative  measurement  of  shock 
sensitivity  (2,42)  the  charge  diameter  effect  on  the  critical  initiating 
shock  was  evident.  This  raises  the  question  of  what  bare  charge  diam¬ 
eter  would  result  in  behavior  equivalent  to  that  exhibited  by  the 
explosive  in  the  regular  gap  test  confinement.  We  showed  some  time 
ago  that  this  confinement  on  cast  Comp  B  gave  a  gap  pressure  predicted 
for  a  bare  charge  of  75  mm  diameter  (43).  The  prediction  could  be 
ma.2  from  a  linear  relationship  established  between  50#  pressure  (P  ) 
and  the  reciprocal  charge  diameter.  The  effective  diameter  of  the 
confined  charge  was  determined  by  replacing  the  confining  material  with 
an  equal  mass  of  explosive,  scill  in  the  cylindrical  configuration. 

If  this  relationship  holds  for  other  explosives  besides  Comp  B,  the 
effective  diameter  dfi  would  probably  have  much  the  same  value  because 
the  high  density  steel  has  much  more  effect  than  the  lower  density 
explosive  in  determining  d  . 

v5 

In  recent  work,  we  have  found  that  a  fine  (ca.  7u)  ammonium  per¬ 
chlorate  dead -presses  in  the  gap  test  at  pQ  *  1.53  g/cc,  i.e.,  the 
material  becomes  subcritlcal  for  propagating  detonation.  Another  fine 
ammonium  perchlorate  (ca.  10p)  exhibited  a  critical  diameter  of  8  mm 
at  pQ  »  1,58  g/cc  (44).  Again  the  gap  confinement  resulted  in  a 
sensitivity  behavior  exhibited  by  a  75-80  mm  unconfined  charge. 

Cast  TNT  lies  between  this  pressed  charge  and  cast  Comp  B  in  its 
shock  sensitivity  behavior.  It  seems  reasonable  to  assume  that  de  of 
cast  TNT  in  the  gap  test  confinement  will  be  about  76.2  mm  (3  inches), 
and  this  value  has  been  used  in  estimating  the  ratio  d  /d  of  Table  1. 

w  w 

Although  the  regular  confinement  results  in  approximately  the  same 
d  for  different  explosives,  the  diameter  effect  on  P  (or  P.)  differs 
markedly  with  the  sensitivity  of  the  explosive.  Thus  -(AP^/Ad  )  is 
35,  l4o,  and  2*10  kbar  cm,  respectively,  for  cast  DINA,  cast  Comp  B, 
and  creamed  cast  TNT.  The  gradient  is  determined  from  P  values 
measured  on  confined  and  unconfined  charges  (d  =  3*81  cm)  with 
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de  —  76  nan  for  the  former.  (See  Table  Al). 

There  are  several  ways  of  estimating  the  initiating  pressures  to 
be  expected  for  confined  charges  of  TNT  Casting  1  and  3  (See  Table  1) 
from  the  measured  values  on  the  unconfined  charges.  The  simplest  is 
to  use  the  linear  relationship  between  gap  thickness  for  the  confined 
and  unconfined  charge.  This  is  indicated  when  the  data  for  cast  DINA, 
Comp  B,  and  TNT  of  Table  Al  are  plotted  in  Fig,  Al.  It  should  be 
mentioned  that  two  pairs  of  pentolite  values  (5,^3)  do  not  fall  on  the 
cur/e,  but  the  range  covered  by  these  two  values  and  six  other  values 
from  confined  samples  is  indicated.  The  upper  part  of  this  area  lies 
quite  close  to  the  curve,  and  makes  it  probable  that  if  the  pentolite 
composition  (PETH/TNT,  50/50)  and  casting  had  been  as  well  controlled 
as  those  of  Comp  B,  its  data  too  would  fall  on  the  curve.  (Many  tests 
were  made  on  pentolite  castings  because  of  difficulties  encountered  in 
trying  to  control  the  slurry  viscosity.) 

The  upper  curve  of  Fig.  Al  shows  the  relation  between  50$  gap 
thickness  of  the  unconfined  (d  *»  3.81  cm)  and  confined  (d  ~  7.6  cm) 
charges.  The  lower  curve  is  drawn  parallel  to  the  upper  through  a 
single  point  in  which  the  unconfined  charge  had  d  =  5.08  cm.  The 
single  point  was  obtained  by  extrapolation  of  the  Comp  B  data  of 
Table  AL.  (For  the  three  measured  values  of  Table  Al,  the  50$  gap 
thickness  varies  linearly  with  the  charge  diameter.) 
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A1  EFFECT  OF  CONFINEMENT  ON  50  %  GAP  IN  LSGT 


NOLTR  70-73 


REFERENCES 


1.  H.  W.  Hubbard  and  M,  H.  Johnson,  J.  Appl.  Phys.  22.,  765  (1959)* 

2.  S.  M.  Brown  and  E.  G.  Whitbread,  Les  Ondes  de  Detonation,  Edition 

du  Centre  National  de  la  Recherche  Scientifique,  Paris,  1962, 

pp  69-80. 

3.  F.  E.  Walker  and  R.  J.  Wasley,  Explosivstoffe  1969,  17(1),  9-13* 
See  also  Univ.  Cal.  Lawrence  Radiation  Laboratory  Preprint  0CRL- 
70891,  26  Jan  1958. 

4.  D.  Price  and  I.  Jaffe,  ARS  Journal  2JLj  595  (1961). 

5.  D.  Price,  I.  Jaffe  and  G.  E.  Roberson,  Ind.  Chim.  Beige  1967. 

32  (Special  No.),  pp  506-510. 

6.  N.  L.  Colebum  and  T.  P.  Liddiard,  Jr.,  J.  Chan.  Phys.  44,  1929 

( 1966) .  — 

7.  I.  Jaffe,  J.  Toscano,  and  D.  Price,  Naval  Ord.  Lab.  Rept. 

NOLTR  64-66  (1964). 

8.  M.  Kamegai  and  J.  0.  Erkman,  "Numerical  Analysis  of  a  Diverging 
Shock  Wave  in  Cylinders  of  Plexiglas",  Preprints  of  5th  Detonation 
Symposium. 

9.  I.  Jaffe  and  A.  R.  Clainnont,  Naval  Ord.  Lab  Rcpt.  NOLTR  65-33 
(1965). 

10.  S.  J.  Jacobs,  T,  P.  Liddiard,  and  B.  E.  Drimmer,  Ninth  Symposium 
(International)  on  Combustion,  pp  517~529,  Academic  Press, 

New  York,  1963. 

11.  D.  Price  and  F.  J.  Petrone,  J.  Appl.  Phys.  35,  710  (1964). 

12.  A.  W.  Campbell,  W.  C.  Davis,  and  J.  R.  Travis,  Phys.  Fluids  4, 

498  (1961). 

13.  G.  K.  Adams,  J.  R.  Holden,  and  E.  G.  Whitbread,  Ind.  Chim.  Beige, 

2  (Suppl. )  pp  216-221  (1959). 

14.  L.  B.  Seely,  J.  G.  Berke,  and  M.  W.  Evans,  AIAA  Journal,  £,  2179 
(1967). 

15.  A.  N.  Dremln,  0.  K.  Rozanov,  S.  D.  Savrov,  and  V.  V.  Yakushev,, 
Combustion,  Explosion  and  Shock  Waves,  2.,  6  (1967). 

16.  D.  Price,  Eleventh  Symposium  (International)  on  Combustion 
pp  693-702,  The  Combustion  Institute,  Pittsburgh,  1967. 


34 


NOLTR  70-73 


17*  D.  Price  and  A.  R.  Clainnont,  Twelfth  Symposium  (International) 
on  Combustion,  pp  761-770,  The  Combustion  Institute,  Pittsburgh, 
1969. 

18.  L.  B.  Seely,  Proc,  Fourth  Electric  Initiator  Symposium  at 
Franklin  Inst.,  Philadelphia,  1963.  Paper  27  of  Rept.  EIS-A2353. 

19.  V.  A.  Vasilyev,  L.  G.  Bolkhovitinov^  and  K.  D.  Khristoforov 
Fizika  goreniya  i  vzryva,  6o4  (1907). 

20.  A.  F.  Belyayev  and  R.  Kh.  Kurbangalina,  Russian  J.  Phys.  Chem. 

2fb  205  (I960). 

21.  A.  N  Dremin  and  V.  S.  Trofimov,  Zhur.  Prikl.  Mek.  i  Tekl.  Fiz. 

No.  1,  126  (1964).  Through  U.  S.  Dept.  Comm.  AD  615211  PP  205-216 
(Mar  1965). 

22«  J.  W.  Enig  and  F.  J.  Petrone,  Phys.  Fluids  398  (1966). 

23.  I.  M.  Voskoboynikov,  V.  M.  Bogomolov,  and  A.  fa.  Apin,  Fiz. 
Goreniya  i  Vzryva  4(1),  45  (1968)  through  JPRS  47007 . 

24.  V.  K.  Bobo lev,  Dokl.  Akad.  Nauk  SSSR  51,  789  (1947).  Translated 
by  U.  S.  Joint  Publications  Res.  Service,  JFRS:402o 

25.  A.  F.  Belyayev  and  M.  K.  Sukoyan,  Combustion,  Explosion,  and 
Shockwaves,  ^(1)  11(1967). 

26.  A.  N.  Dremin,  V.  M.  Zaitzev,  V.  S.  Ilyukin,  and  P.  F.  Pokhil, 
Eighth  Symposium  (International)  on  Combustion  pp  611-619, 
Williams  and  Wilkins,  Baltimore,  1962. 

27.  A.  Ya.  Apin  and  L.  S.  Stessik,  Sbomik  "Fizika  Vzryva"  No.  3, 

1955  USSR  At. ad.  Sci.  Through  Ref  (26). 

28.  L.  N.  Stesik  and  L.  N.  Akimova,  Russ.  J.  Phys.  Chem.  148 
(1959). 

29.  G.  G.  Rempel,  Zh.  Prikl.  Mekh.  Tekh.  Fiz.  1969.  (2),  83-86. 

30.  W.  B.  Cybulski,  W.  Payman,  and  D.  W.  Woodhead,  Proc.  Roy.  Soc. 
(London)  197 A  51(1949). 

31.  I.  Jaffe  and  D.  Price,  ARS  Journal  ^2,  1060  (1962). 

32.  D.  Price  and  A.  R.  Clairmont,  unreported  data. 

33.  M.  J.  Urizar,  E.  James,  Jr.,  and  L.  C.  Smith,  Phys.  Fluids  4, 

262  (1961). 

34.  He  Eyring,  R.  E.  Powell,  G.  H.  Duffy  and  R.  B.  Parlin,  Chem. 

Rev.  4£,  69(1949). 


35 


NOLTR  70-73 


35.  W.  E.  Gordon,  Fourth  Symposium  (International)  on  Detonation 

pp  179-157,  ACR-126,  U.  S.  Gov.  Print.  Office,  Washington,  1967. 

36.  A.  N.  Dremin  and  K.  K.  Shvedov,  Zh.  Prikl.  Mekh.  Tekh.  Fiz.  1964 
(?),  154. 

37*  N.  Griffiths  and  J.  M.  Groocock,  J.  Chem.  Soc.,  I960.  4154. 

38.  H.  Bernier,  J.  M.  Lezaud,  and  F.  Prouteau,  Ind.  Chim,  Beige  1967 , 
32  (Special  No.),  pp  582-586. 

39.  C.  Fauquigon  aod  R.  Cheret,  Twelfth  Symposium  (International)  on 
Combustion,  pp  745-751,  The  Combustion  Institute,  Pittsburgh,  1969 

40.  A.  W.  Campbell,  W.  C.  Davis,  J.  B.  Ramsey,  and  J.  R.  Travis,  Phys. 
Fluids  4,  511  (1961). 

41.  J.  Rcth  and  J.  H.  Blackburn,  "The  Effect  of  Initial  Temperature 
on  the  Shock  Sensitivity  of  Granular  Explosives",  Annual  Report 
of  Stanford  Research  Institute,  28  Aug  1967. 

42.  W.  R.  Marlowe  and  I.  C.  Skidmore,  Proc.  Roy.  Soc.  (London)  246a, 
284  (1958). 

43.  D.  Price  and  I.  Jaffe,  AIAA  Journal  1,  299  (1963). 

44.  D.  Price,  A.  R.  Clairmont,  and  I.  Jaffe,  Combustion  and  Flame  11, 
415  (1967). 


36 


UNCLASSIFIED 


DOCUMENT  CONTROL  DATA  -  R  &  D 

Snofiit  *  Ia\  .it. on  »f  7itfi\  f  .mi  /.mi/.-*  ii.,  •*»  .i-u-f  Lr  *»6<fi  >fir  rrpo  t  r»  (  l,ti.\ihrdl 


i  O  w  _,n**ing  *  C  •  *  *  •  1  *  irp.rMfp  d.if/inri  M  Cgh.N  Cl*r.su  i'  a>'0». 

U.  S.  Naval  Ordnance  Laboratory  UNCIA 

White  Oak 

Sliver  Spring.  Maryland  20910 _ W0WE _ _ _ 

»  o  ^  p  O  p  r  r  •  t ;  j 

Shock  Sensitivity,  A  Property  of  Many  Aspects 


S  TMGfl.ii  *bt r_\t  nmmr,  middle  initial.  Im%l  tuzmr > 

Donna  Price 


t  repor r  date 

15  July  1970 


tm.  CONTRACT  OR  GRANT  NO 

MAT  03L  OOO/fcOll  01  01 

h.  PNOJEC  T  NC 


7 M.  TOTAL  NO  Of  PACES  76.  NO  OF  REES 

33  ^ 


ORIGINATOR'S  RE°OR  T  NUMBER(S) 


N0LTR  70-73 


9t.  OTHER  REPORT  »:OiSJ  Mny  other  numbers  thml  may  be  ms  signed 
this  report  > 


1C  DISTRIBUTION  STATEMENT 

This  document  has  been  approved  for  public  release  and  sale,  its 
distribution  is  unlimited. 

ii  supplementary  notes 

*2  SPONSORING  Ml  LI  T  AR  Y  ACTIVITY 

13  abstrac t 

Naval  Material  Command 

(U)  ‘Shock  sensitivity  of  an  explosive  includes  its  threshold  for 
propagation  of  steady  state  detonation  at  the  critical  diameter  (d  ) 

V 

as  veil  as  its  numerous  thresholds  for  initiation  of  detonation  under 
different  transient  conditions.  Data  for  TNT  charges  of  differing 
degrees  of  homoge  Telty  are  used  to  show  a  continuous  variation  of  d 

C 

with  the  initiating  pressure  P^  measured  with  the  N0L  large  scale  gap 
test.  An  example  of  a  critical  curve  is  constructed  in  the  pressure- 
time  plane;  it  runs  from  the  threshold  conditions  at  d  through  those 
at  P^.  It  is  suggested  that  all  other  initiating  pressures  measured  fc 
the  same  charge  would  also  fall  on  this  curve  which  illustrates  the 
relationship  between  the  different  threshold  values. 
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